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ABSTRACT 

A new scenario for the emission of high-energy gamma-rays from dark matter annihilation around 
massive black holes is presented. A black hole can leave its parent halo, by means of gravitational 
radiation recoil, in a merger event or in the asymmetric collapse of its progenitor star. A recoiled black 
hole which moves on an almost radial orbit outside the virial radius of its central halo, in the cold dark 
matter background, reaches its apapsis in a finite time. Near or at the apapsis passage, a high-density 
wake extending over a large radius of influence, forms around the black hole. Significant gamma-ray 
emission can result from the enhancement of neutralino annihilation in these wakes. At its apapsis 
passage, a black hole produces a flash of high-energy gamma-rays whose duration is determined by the 
mass of the black hole and the redshift at which it is ejected. The ensemble of such black holes in the 
Hubble volume is shown to produce a diffuse high-energy gamma-ray background whose magnitude 
is compared to the diffuse emission from dark matter haloes alone. 
Subject headings: cosmology: dark matter 
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1. INTRODUCTION 

We study the effect on dark matter (DM) distribution 
of black holes (BHs) that are ejected from their parent 
haloes but remain on bound orbits outside the virial ra- 
dius. Black holes can be recoiled due to anisotropic emis- 
sion of gravitational radiation (Bckenstcin 1973, Fitchctt 
1983, Merritt et al 2004). As an ejected black hole (BH) 
moves in the dark matter-dominated universe, a wake 
forms whose density is maximum downstream along the 
symmetry axis. The wake density and its extent, i.e. 
the zone of influence of the BH, increase with decreasing 
velocity of the BH and velocity dispersion of the DM en- 
vironment. A black hole on bound orbit round its parent 
halo comes to rest at apapsis of its highly radial orbit in 
a finite time. In a cold dark matter background with a 
velocity dispersion of the order of cm/s, a BH near its 
apapsis has a large radius of influence. If DM were to 
consist of self-annihilating particles such as neutralinos, 
significant 7-rays would be emitted from these wakes. 
On its apapsis passage, a BH is shown to produce a flash 
of high energy 7-rays whose duration depends on the 
mass of the BH and the redshift at which it is ejected. 
We also estimate the total isotropic diffuse gamma ray 
background emitted from the ensemble of cosmological 
wakes in the Hubble volume, assuming that only a small 
fraction of BHs leave their haloes and compare it to the 
background from DM haloes alone. 

There arc two main assumptions made in this work. 
First, we assume that outside haloes, the cold dark mat- 
ter environment is approximately homogeneous: an as- 
sumption which is more valid at high redshifts. Sec- 
ond, to evaluate the wake density we assume a constant- 
velocity approximation, whose validity is questionable 
and remains to be confirmed by future N-body simula- 
tions. 
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2. DENSITY OF THE WAKE 

The expression for the density enhancement due to a 
moving point mass in a thermal environment is (Danby 
& Camm 1957, Griest 1988) 
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sity of the environment of the BH, A/bh is the mass of 
the BH moving with velocity Vbh and the distance r 
and angle 9 are the radial distance measured from the 
BH position and the angular position is measured from 
the symmetry axis with 9 = upstream in front of the 
BH and 9 = tt downstream on the symmetry axis behind 
the BH. 

When the velocity dispersion of the medium is very 
low in the limit Vbh >> cdm, the integral in ([1]) be- 
comes highly oscillatory and difficult to evaluate. The 
expression for the density enhancement in this limit 
(Vbh >> ctdm) is (see Swcatman & Heggie 2004 for de- 
tails, also Sikivic & Wick 2002 for a different approach) 
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2 vTTcos^ n cosa -I- n^) . the density ([2|) along the sym- 
metry axis (9 = tt) attains the maximum value: p\axis ~ 

^^^nGMBH/r for 9 = n and a^^ « GMbb/t. 
The wake density is independent of the velocity of the 
BH, downstream along the symmetry axis due to the fi- 
nite velocity dispersion of dark matter. 

3. DENSITY ENHANCEMENT AND RADIUS OF 
INFLUENCE: THE ROLE OF BH VELOCITY 

The zone of influence of BH decreases with increas- 
ing velocity and velocity dispersion of its environment as 
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Fig. 1. — The change in the wake overdensity jlj with radial 
distance r from BH in a hot environment (lower curve) and cold 
environment (upper curve). The lower and upper envelopes to the 
black bands correspond to 6 = and d = -k respectively. The solid 
(rod) curves are for stationary BHs ((Sjl. 
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Fig. 2. — The density |(T]| for Vbh < o^DM and l(2} for Vbh >> 
(TDM are shown for a BH moving in a cold (ctdm = 0.03cm/s) [left 
panel] and hot {cr-aM = 2km/s) [right panel] environment. 

shown in Figs. [1] and [2) In Fig. [1] produced by numeri- 
cal integration of equation ([1]), the red middle curve for 
Vbh = shows that this approximation provides a good 
description of the average wake density of slowly-moving 
BHs. In the limit as V ^ Q, the density profile of the 
wake ([T]) reduces to 
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where Erfc is the complementary error function and r, 
is the radius of influence of the BH: r, = Gil/BH/cpjyj. 
We emphasis that ([3]) is the limit 1/ ^ of ([T|), and is 
not a unique density profile for stationary BHs. Here, 
we use ([3]) only as an approximation to ([1]) for slowly- 
moving BHs. Fig. [21 shows the dependence of the den- 
sity enhancement on the BH velocity and DJVI velocity 
dispersion. When the BH is moving fast with respect to 
the background, a significant density enhancement only 
arises in a small zone around the symmetry axis (down- 
stream) of the BH. The density enhancement also de- 



creases with increasing velocity dispersion of dark mat- 
ter environment. The highest density enhancement and 
largest radius of influence are achieved for BHs moving 
slowly in a cold background. 

4. DRAG FORCES ON THE EJECTED BH: TIME TAKEN 
TO REACH THE APAPSIS 

The BH remains bound to its central halo if it is 
ejected with a velocity less than the escape velocity (mea- 
sured from the virial radius). Because it is ejected from 
the centre (and also when with a large velocity), the 
BH is on almost radial orbit. The BH initial veloc- 
ity is set as follows. We assume that the halo mass is 
about 2 X 10'* times the mass of the BH (IVladau & Rees 
2001). Thus, for a BH of mass Mbh, the virial radius of 
the halo, from which it was ejected, can be determined 
using Mhaio = 47r/3Avir(z)p(2:)i?^j,.(z) and noting that 
Avir(z) = (187r2 + 82x - 39x^)/n{z) and x = n{z) - 1 

and n{z) = Qmii + z)^/[nmii + + + + z)^] 

(see Bullock et al 2001 for details). Having evaluated the 
virial radius, we can then evaluate the escape velocity 
from the virial radius of the halo, using T4scapo(-z, -A^) ~ 

vir- 

The ejected BH is slowed down by the gravitational 
pull of its parent halo and also by the dynamical friction 
of dark matter background as 
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where E = —Vj^/2 + G A/haio/^vir is the absolute value 
of the energy with which a bound BH leaves the virial 
radius with velocity Vi. Since dynamical friction plays a 
sub-dominant role in braking the BH, the values of ln(A) 
and the background density p marginally affect the value 
of (dl for a BH in its initial outward journey. We have 
set ln(A) to unity and p to po in order to obtain an 
upper value on the time to apapsis. The left panel of 
Fig. [3l produced from numerical integration of equation 
^ demonstrates that the time elapsed, t^, — ti, since the 
BH was ejected till it reaches the apapsis, is far shorter 
than the Hubble time, to, as long as the BH escapes the 
virial radius with a velocity less than the escape velocity. 
This figure also demonstrates that — ti is independent 
of the BH mass, because higher mass BHs need larger 
velocities to leave their more massive haloes. 

5. GAMMA-RAY FLASHES FROM BHS AROUND APAPSES 

Cold dark matter if composed of neutralinos would self- 
annihilate into secondary products, including energetic 
photons (e.g. see Bertone, Hooper, Silk 2005 for a recent 
review). The absolute luminosity, in units of 7S^* of a 
BH of mass M at rcdshift z is: 
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where is the Schwarzchild radius, r, is the radius of 
influence of the BH, is the neutralino mass (^ 100 
Gev/c^), {(Tv) is the interaction cross-section [which we 
fix at 2 X 10~^^ cm'^/s] and Af.y is the number of photons 
produced per annihilation. Note that the integral ([5]) is 
independent of angle 6 for stationary BHs. 

We had found that the wake density of a slowly mov- 
ing BH is well-approximated by the wake density of a 
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Fig. 3. — Left panel: The ratio of initial BH velocity at the 
virial radius, Vi, to escape velocity, versus ratio of time of ejection 
of BH, ti, to when it reaches apapsis, t*, to the time left from the 
moment of ejection to now {to — tt). Right panel: Minimum time 
the BH spends around the apapsis when its luminosity dominates 
over that of its parent halo [using expression (IIOI I in ^]. The total 
time, tt, taken to reach the apapsis from the virial radius is also 
shown by the dotted horizontal lines for different redshifts, for BHs 
ejected with near escape velocity. 

stationary BH. By inserting ^ [keeping only the first 
term] in ^ we obtain the analytic expression (in units 

of 7 s~^) 



L 



BH 



1.4 X 10^^(l + z) 



Ah 



BH 



Mrr 



'cutoff 



(6) 



for the absolute luminosity of a BH, where i?cutoff is in 
unit of parsecs. In ([5]) i?cutoff is the radius within which 
the luminosity of the BH is evaluated. 

Next, the BH luminosity is compared to the abso- 
lute luminosity of its central dark matter halo of mass 
2 X IO^Mbh, assuming it has a NFW density profile 
(Navarro, Frenk & White 1997) 
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where n{z) = f7,„(l + zf/in^il + z)^ + + f7fc(l + 
z)2], Sc = (Avir/3)c3/(ln(l + c) -c/(l + c)), and c is the 
concentration parameter, for which we use the fit [10/(1+ 
z)](Mvir/M)"° " (agrees weh with Bullock et al 2001, 
Hennawi et al 2007) and Afhaio = Mvh- is the halo mass 
within the virial radius. The absolute luminosity L ([5]), 
in unit of 7 s^^, of a NFW halo of mass Mhaio can then be 
evaluated using ([7]) in ([5]) [after setting M = Mhaio,''. = 
-Rvir, fs = 0] and fitted by the following functions 
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We can find an analytic expression for i?cutoff by using 
the dynamical friction formula ([4]) and assuming that 



near the apapsis dynamical friction dominates over the 
gravitational pull of the halo. This yields 
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where At* and i?cutoff are defined as the time-duration 
and distance-parcours during which the velocity of BH 
reduces from a fraction a of the background velocity dis- 
persion to zero (at the apapsis), i.e. < Vbh < oia{z) 
where a{z) := (1 + z)aQ and (Tq is present velocity disper- 
sion of DM (approximately 0.03 cm/s for neutralinos). 
The lower-bound on i?cutoff and At* is obtained by re- 
quiring that Lbh/-^nfw > 1 where Lbh is given by ([6]) 
and Lnfw is given by ([8]). These lower bounds are shown 
in Fig.[4]and Fig. [3] respectively. Fig|4]shows that at high 
redshifts most ejected BHs are more luminous than their 
parent haloes whereas at low redshifts a much higher 
^cutoff is required. Similarly, we put a lower limit on the 
parameters a 



a > 



7 X 10^ /M, 
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for z > 1 and for z < 1 we can use the very similar 
expression a > (10^/(1 + z)3/'^)(Mbh/M0)^/'' which are 
shown in Fig. H) The two panels of this figure show that 
less massive BHs at higher redshifts are the most lumi- 
nous. 

6. DIFFUSE GAMMA-RAY BACKGROUND 

In a cosmological scenario, ejected BHs near their 
apapses pasages, especially those at high redshifts, where 
the merger rate is higher, can yield an observable diffused 
background flux. The total flux is given by the integral 

^^i^N{M,z)dMdV(z) (11) 
47rr(z)^ 
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where M can be either the BH mass (Mbh) or the halo 
mass (Mhaio), r{z) = Rh ^1 — y^j) with Hubble ra- 
dius Rh — 4000Mpc, and N{M, z) is the number density 
of the BHs [or haloes in the calculation for NFW haloes] 
and the luminosity of a single BH i(M, z) (or the parent 
NFW haloes at z) is given by ^ [or ([H])] and the volume 
element is dV = simp r(z)'^dr{z) dtp dtp. 

The physical number density of the BHs is assumed to 
follow the Press-Schechter formalism (Press & Schechter 
1974, Bower 1991), multiplied by the relative time a BH 
spends at apapsis, i.e. (At*)/to where to is the age of 
the Universe. 

The time spent at apapsis is itself a function of BH 
mass as given by ^ and shown in the right panel of 
Fig. [3] for minimum value of a [obtained by substitut- 
ing (Uni) in (O]. Thus only the fraction At* /to of the 
ejected BHs can be considered to be actually on their 
apapses passage in a Hubble time. We can evaluate the 
total flux by performing the integrals in (jll[) [multiplied 
by Min[At*]/to for BHs]. For spectral index n = — 1 
and Af* ~ 10^^ and for the ensemble of BHs at 
their apapses passages and their central haloes (assuming 
NFW profiles) we obtain <i>NFW ~ 10~^ 7 cm~2sr~^ 
The flux from the BHs is lower than this value (since we 
have evaluated our parameters i?cutoff and a by requir- 
ing ivNFW = ^bh) because the BHs only spend a fraction 
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Fig. 4. — left panel: The minimum radius i?cutoff apapsis 
required for the BH luminosity to dominate over that of its parent 
halo is shown by the lower solid line. The minimum virial ra- 
dius, iJviri of the halo (corresponding to ejection of a BH of mass 
IOOMq) and similarly the minimum distance from the virial radius 
of the apapsis Rt are also plotted, right panel: the minimum value 
of the parameter o as a function of redshift and mass of the ejected 
BH [expression lllUl l] are shown. 

of time at the apapsis [Min(AtH,)] which yields approxi- 
mately $NFW ~ 10~^^ 7 cm~^sr~^ The flux would 
be attenuated due to interaction of photons which how- 
ever would affect approximately equally <&bh and <&nfw- 
We have assumed that only BHs produced in 3a peaks of 
the density perturbation can undergo effective mergers. 
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We have assumed that all ejected BHs orbit their central 
haloes outside the virial radius; however, were this not 
the case, we do not expect any significant overall decrease 
in the flux which is already underestimated by our mod- 
erate choices of parameters and also by assuming that 
there is only one apapsis passage for a BH. We have also 
ignored the effect of multiple density-enhancement for a 
BH which is on its inward journey through an already 
high-density wake. 

In conclusion, BHs on bound orbits around haloes can 
be powerful sources of high-energy 7-rays, both indi- 
vidually as resolved sources and collectively as diffuse 
background. The results here indicate that the globu- 
lar clusters in the outskirt of our halo or field galaxies 
in our local Universe devoid of central BHs can have or- 
biting BHs which during their apapsis passages would 
produce flashes of high-energy 7-rays, although this ef- 
fect is expected to be most significant at high redshifts. 
The validity of dynamical friction formulae has been very 
rarely studied for radial orbits (Gualandris & Merritt 
2007). The fact that there is no mass- loss makes BHs 
a rare case for dynamical friction theory. Throughout 
this work we have assumed a homogeneous background 
and a constant- velocity approximation, both of these as- 
sumptions are questionable for the problem considered 
here. It thus remains for forth-coming numerical works 
to check the extent of the validity of the present results. 
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